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The quantum analog of Carnot cycles in few-particle systems consists of two quantum adiabatic 
steps and two isothermal steps. This construction is formally justified by use of a minimum work 
principle. It is then shown, without relying on any microscopic interpretations of work or heat, 
that the heat-to-work efficiency of the quantum Carnot cycle thus constructed may be further 
optimized, provided that two conditions regarding the expectation value of some generalized force 
operators evaluated at equilibrium states are satisfied. In general the optimized efficiency is system- 
specific, lower than the Carnot efficiency, and dependent upon both temperatures of the cold and 
hot reservoirs. Simple computational examples are used to illustrate our theory. The results should 
be an important guide towards the design of favorable working conditions of a realistic quantum 
heat engine. 


Introduction - The big energy challenge of this century 
calls for diversified energy research, including a bottom- 
up approach towards energy efficiency. Apart from two 
stimulating implementations of microscale heat engines 
mm, some theoretical aspects as well as possible real¬ 
izations of nanoscale heat engines [MS] have been stud¬ 
ied. For purely quantum heat engines at the nanoscale 
where the working medium may consist of few particles 
only (e.g., few trapped ions 0)> both quantum fluctua¬ 
tions and thermal fluctuations become significant. Gen¬ 
eral understanding of the design of such energy devices 
are also of fundamental interest to nanoscale thermody¬ 
namics [171420] . In particular, as the size of the working 
medium shrinks to a quantum level, one must reexamine 
the implications of the second law of thermodynamics for 
the efficiency of quantum heat engines. To that end, we 
construct and look into the quantum analog of Carnot 
cycles [2ll|22]. 

The construction of the quantum analog of a Carnot 
cycle is not as straightforward as it sounds. Consider 
first the two quasi-static isothermal steps during which 
the working medium is in thermal equilibrium with a 
reservoir. Regardless of the size of the quantum medium, 
its thermodynamic properties can therefore be well de¬ 
fined in the standard sense. As such isothermal steps can 
be directly carried over to the quantum case. However, 
translating the two adiabatic steps of a Carnot cycle into 
a quantum analog is by no means obvious. One intu¬ 
ition [ai m no] is to replace quasi-static adiabatic steps 
in thermodynamics (without heat exchange) by quantum 
adiabatic processes (as defined in the celebrated quantum 
adiabatic theorem [23]). The starting point of this work 
is to formally justify such an intuitive construction by re¬ 
vealing a fundamental reason related to energy efficiency. 

Below we simply call the quantum anolog of a classi¬ 
cal Carnot cycle (two isothermal steps and two quantum 
adiabatic processes) as a quantum Carnot cycle. It is yet 
fundamentally different from a conventional Carnot cy¬ 
cle. During the two adiabatic steps, the working medium 


implementing the quantum Carnot cycle is generically 
not at equilibrium conditions, except for the case in which 
all the energy levels of the working medium are scaled 
by a common factor as a system parameter varies (to be 
elaborated below). Thus, it becomes important to lay out 
general designing principles concerning how the efficiency 
of a quantum Carnot cycle can be optimized, preferably 
using standard definitions of work and heat. The explicit 
optimization conditions are presented below. Our theory 
also shows that in general the optimized efficiency at¬ 
tained by a quantum Carnot cycle is (i) lower than the 
standard Carnot efficiency, (ii) not a simple function of 
TjTy but a function of both T c and Ty, the temperatures 
of cold and hot reservoirs, and (iii) depends on the de¬ 
tailed spectrum of the working medium. These features 
will guide us in the design of favorable working conditions 
of a realistic quantum heat engine. Simple computational 
examples are used to illustrate our theory. Throughout 
this work we do not use recent microscopic interpreta¬ 
tions or definitions of work or heat proposed for quantum 
systems, such as those introduced in Refs. in uni nn [24]. 
Instead, we only assume that heat exchange is zero if the 
working medium is thermally isolated and work is zero if 
the system parameters of the working medium are fixed. 

Efficiency of quantum heat engine cycles and the sec¬ 
ond law - We start with general considerations of a quan¬ 
tum heat engine cycle consisting of two isothermal steps 
and two thermally isolated processes. Figure [l] schemati¬ 
cally depicts such a cycle. There A B and C —> D rep¬ 
resent two isothermal processes during which the quan¬ 
tum medium is always at equilibrium with a reservoir, A 
is assumed to be the only system parameter tunable in 
a cycle opeation, (. E ) is the mean energy of the system, 
B —>• C' and D —>• A' represent two thermally isolated 
and hence unitary processes. The symbols A' and C’ are 
to indicate that right after a unitary process, the quan¬ 
tum medium is in general not at thermal equilibrium. 
States of A! and C' will reach thermal equilibrium states 
A and C after relaxation with a reservoir under fixed 
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FIG. 1. (Color online) A quantum heat enegine cycle consist¬ 
ing of two isothermal steps and two thermally isolated and 
hence unitary steps. A, B , (7, and D represent four equilib¬ 
rium states, C' and A' represent two non-equilibrium states 
at the end of two unitary steps, approaching respectively to 
equilibrium states C and A after a relaxation step initiated by 
contact with the reservoirs. As shown in the text, for a quan¬ 
tum analog of the Carnot cycle, the two thermally isolated 
steps should be two quantum adiabatic processes. 

values of A. 

When the system at the nonequilibrium state C' starts 
heat exchange with the cold reservoir under fixed A = 
Ac, no work is done. Hence (Ec f ) — {Ec) is simply the 
heat dumped (which could be negative). This thermal 
relaxation process is followed by the isothermal process 
from C to D. The total heat Q 0 ut dumped to the cold 
reservoir (Q out > 0 indicates heat flowing from system to 
cold reservoir) is hence contributed by two terms, with 

Qout = T c (Sc ~ S D ) + (Ec) - (Ec), (1) 

where Sc and Sd describe entropy of equilibrium states 
C and D. In the same fashion, the total heat absorbed 
from the hot reservoir, denoted by Q in (Qin > 0 indicates 
heat flowing from hot reservoir to system), is given by 

Qin = T h (Sb - Sa) + (Ea) - {Ea /). (2) 

The efficiency of such a general quantum engine cycle is 
therefore r] q = 1 - Q 0 ut/Qin, i-e., 

_ 1 _ Use ~ S D ) + (Ec) z (E c ) , . 
Vq T h (S B -S A ) + (E A )-{E A .y 

To compare the above efficiency rj q with the Carnot ef¬ 
ficiency rj c = 1 — T c /Th, we first define AS^f^ and 
A5'^ t ^| 4 , namely the total entropy increase of the uni¬ 
verse for the overall process D A and B — C . It is 


straightforward to obtain 

A«c = (S C - S B ) - h(E c ) - (E c >)]-, 

-L c 

A S%? A = (S A - S D ) - ±-[(E A ) - (E A ,)). (4) 

By the second law of thermodynamics, both Aand 
Acannot be negative. Let us now rewrite Eq. ^ 
as: 

„ = _ T c (S B -S D )+T c ASg? c 
Vq T h (S B - S D ) - T h AS %*? A ' 1 1 

Evidently then, if A S)j^ l A and A S^^ l c in Eq. ^ is zero, 
then rj q would reduce exactly to the Carnot efficiency r\ c . 
In general, r] q in Eq. ([5| is seen to be lower than r] c . In 
short, the second law of thermodynamics implies that, 
the efficiency of a quantum heat engine cycle described 
above should be in general lower than, and can only reach 
in exceptional cases, the Carnot efficiency. 

Constructing a quantum Carnot cycle - To construct a 
quantum Carnot cycle, one must specify the two uni¬ 
tary processes B —>> C' and D A !. Reference [3] 

first proposed to consider quantum adiabatic processes 
for this purpose, mainly based on reversibility considera¬ 
tions [25] . Here we show that this intuitive construction 
is correct for a more fundamental reason related to the 
heat-to-work efficiency. 

Before proceeding, we emphasize that adiabaticity in 
a quantum unitary process does not have the key fea¬ 
ture of a thermal quasi-static adiabatic process in the 
Carnot cycle, i.e., the former does not result in equilib¬ 
rium states in general but the latter does. As a result, 
quantities such as temperature and thermodynamic are 
usually ill-defined for states C’ and A'. Consider then 
the expression of rj q in Eq. (|3|. With four equilibrium 
states A, B, C and D specified as in Fig. 0 only {Ec) 
and (. Ea f ) may be varied by choosing different types of 
unitary processes B C' and D A'. For thermally 
isolated processes, there is no heat exchange and as such, 
we have 

{Ec) = {E b ) + {W) b ^ C '-, 

(E a ,) = {E d ) + (W) D ^A’ (6) 

where (W)b->C' an d (W)d^A' represent the average 
work associated with B -A C’ and D -A A !. Remark¬ 
ably, the minimal work principle |26j then takes us to a 
definite choice. In particular, for a quantum state ini¬ 
tially prepared as a Gibbs equilibrium distribution (this 
specific requirement can be loosened) and for fixed initial 
and final A values, a quantum adiabatic process (if im- 
plementable) is the one with the minimal average work. 
So if D ^ A! and B —>► C’ are indeed quantum adia¬ 
batic processes, the minimal work principle ensures that 
the final mean energies {Ec) or {Ea>) are minimized 
for fixed states B and D. Returning to the expression 
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of r) q in Eq. minimized (Ec') and (Ea') then yield 
the highest possible efficiency r\ q . It is for this efficiency 
consideration that the quantum analog of Carnot heat 
engines must consist of two quantum adiabatic steps in 
addition to two isothermal steps. To our knowledge, this 
is an important and previously unknown insight [25] . 

Optimizing efficiency of quantum Carnot cycles - With 
quantum Carnot cycles constructed and justified as 
above, we next seek specific design principles to further 
optimize rj q . The Hamiltonian of the working medium 
(when thermally isolated) is assumed to be H( A) with 
energy levels E n ( A). The values of A at B and D, namely, 
A b and A£>, are assumed to be given. The focus question 
of this study is to show how to choose A at states A and 
C, namely, A^ and A c, such that r] q may be optimized. 

It is interesting to first illustrate this optimization issue 
in systems possessing scale invariance [27l 128] with A. In 
such an exceptional case, 


[En( Ai) - E m (Xi)] = S(X u X 2 )[E n (X 2 ) - E m ( A 2 )]. (7) 

Examples of this situation include a harmonic oscillator, 
with A being as the harmonic frequency, a particle in 
a infinitely deep square-well potential [3], where A can 
be the width of the potential well, or simply a two-level 
quantum system [9]. Consider then the adiabatic step 
from B to C'. The initial state populations are given by 
Pb{ n) = e -^hE n (\ B )/Z b (throughout Z represents equi¬ 
librium partition functions and (3 represents the inverse 
temperature). Upon reaching C", the final populations 
are still given by Pc(n) = e~^ hEn ^ XB ^ /Z B due to the 
assumed quantum adiabaticity. Now given the assumed 
scale-invariance in Eq. 0 , one can always define an ef¬ 
fective temperature T e ff to reinterpret Pc (n ), namely, 

Pc (■ n ) = e“^ E " (AB) /Z B = e“ Aff£ " (Ac) /Z c , (8) 


following: 


Qout = (Ec) - (Ed) ~ k B T c In f^, 

Zc 

Qin = ( E b ) - (Ea'} + k B Th In 

Za 


(9) 


Interestingly, with states B and D fixed, only A^ may 
affect Qout via (Ec) and Zc ; while only A^ may affect 
Qin via (Ea') and Za • That is, to optimize 77 ^, mini¬ 
mizing Q ou t and maximizing Q in can be executed sepa¬ 
rately, which is a considerable reduction of our optimiza¬ 
tion task. For this reason, below we focus on minimizing 
Qout and the parallel result concerning Q in directly fol¬ 
lows. 

Accounting for quantum adiabaticity that maintains 
populations on each quantum level, one has 


(Ec) = e- phE ^ XB) E n (\ c ). ( 10 ) 


Note again that the level populations A-e~^ hEn ^ XB ^ used 
above are in general not an equilibrium Gibbs distribu¬ 
tion associated with H(Xc )• Using Eqs. ^ and ( 10 ), we 
arrive at 


<9Q out 

dX c 


d(Ec) 


+ k B T c 


dQnZc) 


E 


d\c ~ L dX c 

e —f3hE n (X B ) e —P c En{X c ) 


Z B 


Zc 


d_Eu(Xc) 

dX c 

( 11 ) 


The minimization of Q 0u t requires the condition 
OQovt/OXc = 0 , which indicates that 


E 


"g— PhE n (\ B ) 

~Z~b 


e -/3 c E n (A c )- 


dE n (X c ) 

dX c 


( 12 ) 


where /3 e q = 1 /(fc^Teg) = S(X B , Xc)fih- That is, state 
C’ has no difference from an equilibrium state with tem¬ 
perature T e ff and Hamiltonian H(Xc ). If we now choose 
A c to guarantee that T e ^ = T c , then state C’ is already in 
thermal equilibrium with the cold reservoir at T c . The re¬ 
laxation process from C' to C as illustrated in Fig. |T]is no 
longer needed, resulting in Sc = (Ec) = (Ec) and 
hence ASffffc = 0. Exactly the same analysis applies to 
the adiabatic process D —)> A !. That is, by choosing an 
appropriate value of A^, we can set A Sq^ 1 a = 0 . Ac¬ 
cording to the expression of g q in Eq. ([5|, r] q then yields 
the standard Carnot efficiency. This result also offers 
a clear perspective to explain why the Carnot efficiency 
can be obtained in some early studies of quantum heat 
engines [3JCE31. 

We next lift the above scale-invariance assumption and 
proceed with optimizing 77 ^, by optimizing Q ou t and Qi n . 
We first rewrite Q ou t and Q in in Eqs. 0 and 0 as the 


Viewing the linear response in energy to a variation in A 
as a generalized force, we define a general force operator 
j-\ = — • Then the condition in Eq. |l2| can be 

cast in the following compact form, 

(E\ c )c = (Ub^ c P\ c U b ^c)b, (13) 

where Ub^c is the unitary transformation that trans¬ 

forms an arbitrary nth eigenstate of H(Xb) to the Trth 
eigenstate of H(Xc ). That is, the expectation value of a 
generalized force operator at A c over equilibrium state C 
should be identical with that of a mapped force operator 
over equilibrium state B. Needless to say, the condition 
for Q in to be maximized is given by 

(E\ a )a = (Ud^ A E\ a Ud^a)d, (14) 

where Ud^a transforms an arbitrary 72 th eigenstate of 
H(Xd) to the Trth eigenstate of H(Xa )• 
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Unlike a previous interesting suggestion [29] . the two 
explicit conditions in Eq. (13) and (14) to optimize rj q 
are not about matching the mean energy between states 
C' and C ( A' and A). Attempts to match information 
entropy between states C' and C ( A' and A) do not op¬ 
timize Tj q , either. Rather, the conditions found here are 
about a more subtle and more involving matching of the 
expectation values of some generalized force operators 
through equilibrium states. We now take Eq. (12) as 
the example to digest the optimization conditions. For 
the exceptional case of a scale-invariant medium, due 
to the existence of a /3 e ^ = /3 C at C", one can achieve 

0 for arbitrary n. Then 


0 -(3 h E n (\ B ) 


~/3 c E n (\ c ) 




Eq. (112]) can be easily satisfied, independent of the details 
of dK Qx Xc ^ • For a general working medium, the condition 
of Eq. (12) may be still satisfied after setting the sum of 

all the terms dE Q^ c ^ e /3 ^ (Ag) _ e Ac) zero . 

Numerical examples - We adopt a simple model sys¬ 
tem that is not scale-variant with A, with E n ( A) = 
An+cm 2 +const (all variables in dimensionless units). For 
our purpose here there is no need to specify the explicit 
form of the Hamiltonian. If a is comparable to A, then 
the ratio [Fk(Ai) - E rn (X 1 )]/[E n (X 2 ) - E rn (X 2 )} does de¬ 
pend strongly on n and m, a clear sign of breaking the 
scale invariance. Cases of a very small a would resemble 
the behavior of a harmonic oscillator at low tempera¬ 
tures. To guarantee quantum adiabaticity, we exclude 
cases with level crossings. This is achieved by requiring 
A > — a such that E n > E m if n > m. Other physical 
considerations for the cycle to operate as a meaningful 
heat engine suggests that A a should be the largest and 
A c the smallest among Aa, Ad, Ac, and Ad- Our com¬ 
putational details confirm that minimization of Qout and 
maximization of Q in indeed occur precisely at those loca¬ 
tions predicted by Eq. (13) and (14). Optimization of r\ q 
as we outlined above theoretically is hence indeed doable. 

The rather specific conditions to optimize rj q (under 
fixed A b and A d) indicates that the r] q thus optimized 
will be highly system specific. To see this, we present in 
Fig. [^optimized rj q as a function of Ad (Ad) with fixed 
Ad (Ad), under different temperatures T c and 7k. From 
Fig. |2(a)| it is seen that the optimized r} q can be way 
below, but nevertheless quickly approaches, the Carnot 
efficiency r) c as Ad increases. This is because a larger 
Ad leads to an even larger A a due to the optimization 
requirement, both facts pushing the system closer to a 
scale-invariant system under fixed temperatures. Note 


also that even though the three r] q curves in Fig. 2(a) 


are for the same ratio 7k/7k, their rj q values are much 
different. This shows that the optimized r\ q is no longer 
a simple function of T c /7k, but a function of both T c and 
Th. Figure 2(b) shows that our optimized r\ q may not 
be always a monotonous function of Ad with a fixed Ad • 
Interesting effects of T c and 7k under a common 7k/7k 
are again observed there. The shown ranges of Ad or Ad 


in Fig. [2] vary with the chosen temperatures because we 
exclude level crossing situations. The sharp change of r] q 
in Fig. |2(a)| [Fig, 2(b)| with a decreasing (increasing) of 
A d (Ad) under a given A d (Ad) is simply because the 
optimized cycle is about to cease to operate as a heat 
engine (which requires Q[ n ^ Qout ^ o)* 



(a) 



(b) 


FIG. 2. (Color online) Behavior of heat-to-work efficiency rj q 
optimized under given Ad, Ad, T c and Th , with A a and Ac 
chosen according to Eq. (13) and Eq. (a) r\ q vs Xd if 

A b = 6.0 for three sets of T c and Th. and (b) rj q vs A b if 
A d — 6.0, also for three sets of T c and Th. In all the shown 
cases the Carnot efficiency rj c *= 0.5 (top dashed curve). All 
variables are in dimensionless units, and the energy levels of 
the working medium are assumed to be E n ( A) = An + an 2 
with a — 0.1. 


Discussions and conclusions - Several previous studies 
investigated quantum Otto cycles 

[iHUCElEn] consist¬ 
ing of two thermally isolated steps and two isochoric pro¬ 
cesses that are simply relaxation processes with a hot or 
a cold reservoir. It is clear that such quantum Otto cycles 
can be regarded as a special case of the quantum heat en¬ 
gine cycles considered here, without the isothermal pro¬ 
cess A —>> B or C -A D. That is, by setting A a = Ad and 
Ac = Ad in Fig. 0 we obtain the quantum Otto cycles. 
One can now also justify the use of quantum adiabatic 
steps to construct energy efficient quantum Otto cycles, 
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using the minimum work principle again [26]. But more 
importantly, because in our efficiency optimization under 
fixed A b and A£>, the obtained A a (Ac) in general differs 
from A b (Ad), one deduces that the optimized r] q here is 
in general higher than the efficiency of the corresponding 
quantum Otto cycles. This fact strengthens the impor¬ 
tance of quantum Carnot cycles we have justified and 
optimized. 

Our analysis of the two adiabatic steps does not re¬ 
ally demand the two steps to be executed slowly. That 
is, so long as the final-state populations are consistent 
with those expected from the quantum adiabatic theo¬ 
rem, then all our results will be valid. This understanding 
encourages the use of shortcuts to adiabaticity [28j 30 
[33] or even an optimal control approach |34] to implement 
the quantum Carnot cycles within a shorter time scale, 
thus boosting the heat engine power. 

In conclusion, using minimal assumptions about the 
concepts of work and heat in few-body systems, we have 
shown how to construct and optimize the quantum ana¬ 
log of Carnot cycles at the nanoscale. The heat-to-work 
efficiency can be optimized if two conditions regarding 
some generalized force operators evaluated at some equi¬ 
librium states are met. In general the optimized efficiency 
is system specific, lower than the Carnot efficiency, and 
dependent upon both temperatures of the cold and hot 
reservoirs. 
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